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ABSTRACT 

Two  basic  approaches  to  the  problem  o£  estimating  bearing  in  the  ocean  are 
compared.  These  approaches,  both  of  which  utilize  signal  phase,  rather  than  ampli¬ 
tude,  information,  are  the  traditional  closed-loop  method  and  a  newer  open-loop 
concept.  The  relative  theoretical  performance  is  determined  in  terms  of  the  stan¬ 
dard  deviation  of  fluctuation  in  the  bearing  estimate,  normalized  for  universal¬ 
ity.  Bias  considerations  are  neglected  for  simplicity.  When  the  resulting  closed- 
and  open-loop  performance  curves  are  compared,  it  is  found  that  there  is  little 
difference  between  them  at  high  signal-to-noise  ratios. 

K 


HIGHLIGHTS  OF  REPORT 


■  Bearing  information  imbedded  in  differential  signal  phase  is  identified. 

■  Basic  configurations  of  closed-loop  and  open- loop  bearing  estimators  used  in 
the  study  are  defined. 

■  The  phase  detection  functions,  although  different  in  appearance,  prove  to  be 
mathematically  identical. 

■  Phase  statistics  are  derived. 

■  The  location  and  nature  of  the  error  differ  significantly  between  the  two 
estimators. 

■  Performance  of  the  two  estimators  is  evaluated  and  compared. 

■  The  theoretical  performance  of  the  two  estimators  is  the  same  at  high 
signal-to-noise  ratios. 

■  The  decision  as  to  which  estimator  to  choose  should  be  based  on  considera¬ 
tions  other  than  theoretical  performance. 
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RECOMMENDATIONS  FOR  FUTURE  STUDY 

■  This  study  has  made  sevecal  simplifying  assumptions: 

•  No  source  motion 

•  Source  exactly  on  beam  main  response  axis 

•  Noise  in  the  two  half  beams  statistically  independent  of  each  other. 

A  more  realistic  comparison  would  result  if  these  assumptions  were  to  be 
relaxed.  This  is  recommended  for  a  future  study  task  or  combination  of 
tasks. 

■  This  study  has  examined  the  use  of  signal  phase  to  estimate  bearing.  It  is 
recommended  that  a  similar  study  be  initiated  to  examine  the  use  of  signal 
amplitude  instead:  fundamentally,  the  performance  of  a  bearing  interpola¬ 
tion  scheme.  Other  studies  of  this  type  have  been  performed,  but  not  neces¬ 
sarily  asking  the  questions  that  need  to  be  asked  for  this  application. 


v 
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1.0  IWHOOqCTIOM 


This  is  a  technical  ceport  prepared  under  the  "DLT  Anomaly  Study"  performed  by 
the  author  under  the  technical  direction  of  Peter  J.  Barnikel,  NUSC/NLL,  under  con¬ 
tract  N00140-78-6021-0052. 


The  question  addressed  in  the  present  report  is  the  following:  Is  there  a 
fundamental  difference  in  theoretical  performance  between  a  closed-loop  bearing 
estimator,  typified  by  the  DIFAR-Like  Tracker  (DLT) ^ ,  on  the  one  hand,  and  an 
open-loop  bearing  estimator  on  the  other  hand,  both  estimators  making  use  of  signal 
phase,  rather  than  amplitude,  information? 

It  should  be  mentioned  in  passing  at  this  point  that  studies2-4  have  been 
performed  on  the  subject  of  bearing  estimation  (interpolation)  using  amplitude 
information.  Despite  this,  it  would  be  wise  to  perform  such  a  study  in  the  future, 
asking  and  answering  the  right  questions  pertinent  to  the  application  at  hand.  In 
particular,  performance  should  be  compared  with  the  performance  results  obtained  in 
the  present  report. 

This  report  is  quite  theoretical,  making  many  simplifying  assumptions  and 
working  with  normalized  quantities  to  arrive  at  a  performance  compar ison  between 
idealized  classes  of  systems.  At  no  time  is  any  information  on  the  absolute  per¬ 
formance  of  any  real  system  derived  or  provided;  that  would  be  beyond  the  scope  of 
the  present  study. 

Since  the  bearing  estimators  under  consideration  use  phase  information,  the 
second  section  of  this  report  is  devoted  to  an  examination  of  the  information  car¬ 
ried  in  signal  phase.  The  third  section  shows  the  basic  configurations  of  the 
closed-loop  and  open-loop  bearing  estimators  under  study  in  this  report.  The 
fourth  section  analyzes  the  respective  phase  detection  implementations.  The  fifth 
section  presents  the  phase  statistics  derived  in  Appendix  A.  The  sixth  section 
examines  the  location  and  nature  of  the  error  in  the  two  estimators.  Finally,  the 
performance  of  the  estimators  is  evaluated  and  compared  in  the  seventh  section. 


1.  References  are  listed  in  Section  8 
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2.0  MASK  IMTO8NATI0B 


The  two  bearing  estimation  concepts  being  compared  in  this  report  both  make 
use  of  information  contained  in  the  phase  of  the  arriving  signal.  In  this  section, 
the  information  content  of  the  signal  phase  is  examined. 

First,  the  signal  is  described!  It  will  be  assumed  in  this  report  that  the 
signal  waveform  is  a  pure  sinusoid  and  that  the  signal  wavefront  is  plane.  That 
is,  the  radiating  source  is  assumed  to  be  in  the  far  field  and  the  propagation  is 
assumed  to  be  ideal,  introducing  no  multipath  structure  and  no  fluctuation,  either 
in  amplitude  or  in  phase.  In  general,  relative  motion  can  exist  and  have  an 
effect)  in  the  end,  however,  relative  motion  will  be  assumed  not  to  exist. 

Second,  the  receiving  array  is  described:  Por  both  bearing  estimation  con¬ 
cepts,  an  array  of  hydrophones  split  into  left  and  right  halves  is  used.  The  base¬ 
line  distance  between  phase  centers  of  these  halves  is  called  the  effective  dipole 
spacing  (d) . 

The  signals  appearing  in  the  left-  and  right-half  beam  outputs  can  be  written 
as  follows: 


SL(t)  »  A  cos 
SR(t)  *  A  cos 


[2TTf0t  +*L(t)] 
[2nf0t  +  *R(t)] 


(2-1) 

(2-2) 


where 


A  *  constant  signal  amplitude,  usually  the  same  in  both  half  beams 
(although  not  a  requirement) , 

£q  ■  source  frequency  of  signal, 

$L,R  ■  phase  terms  in  left  and  right  half  beams  respectively,  further 
described  as  follows: 
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where  k 


« 


L,R 


j-ritt;  +  -|-  £sin  e1(t)  -  sin  0Q  J  | 


acoustic  wave  number: 


where 


(2-3) 


(2-4) 


X.  *  acoustic  wavelength 

c  *  speed  of  sound 

Continuing, 

*  range  to  source 

■  bearing  to  source 

0Q  *  bearing  direction  of  array  main  response  axis  (MRA) 

■  constant  additive  phase  shift. 

Now  the  information  contained  in  the  phase  terms  given  by  (2-3)  will  be  dis¬ 
cussed.  To  start,  the  constant  additive  phase  shift  (4>Q) ,  of  course,  conveys  no 
information  at  all. 

Next,  source  motion  needs  to  be  discussed.  The  DLT  study*  analyzes  a 
closed-loop  bearing  estimator;  in  that  analysis,  source  motion  is  taken  into 
account.  One  practical  effect  of  substantial  source  motion  is  the  creation  of  a 
need  for  automatic  frequency  tracking  (AFT).  The  range  term  (r^)  provides  the 
information  required  for  the  AFT  function:  Doppler  shift  is  proportional  to  the 
time  derivative  of  r^(t).  The  present  study  is  cast  into  a  different  frame  of 
reference:  Two  bearing  estimators,  one  closed-loop  and  the  other  open-loop,  are 
being  compared.  To  simplify  the  comparison,  source  motion  is  assumed  not  to 
exist.  Therefore,  in  this  case,  the  range  term  (ri)  is  constant  and  conveys  no 
more  information  than  does.  Consequently,  the  AFT  function  is  no  longer  nec¬ 
essary,  hence  will  receive  no  further  consideration  in  this  report.  Furthermore, 
the  input  bearing  (0j)  is  also  constant. 
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So  focus  on  the  informative  ton*  and  to  eliminate  useless 
called  to  the  differential  phase i 

'S  s  #1,  -  $  £ 

«  -  Kd  (sin  Oj  -  sin  0O) 


terms,  attention  is 


(2-5) 

(2-6) 


This  last  equation  requires  further  examination.  For  a  closed-loop  bearing  estima¬ 
tor  the  bearing  error 


©e  *  ©i  -e, 


(2-7) 


will  be  relatively  small  if  the  input  signal-to-noise  ratio  (SNR)  is  not  low,  par¬ 
ticularly  when  source  motion  is  absent.  For  the  open-loop  bearing  estimator,  this 
statement  is  not  necessarily  true  (this  matter  will  be  discussed  at  greeter  length 
in  a  later  section),  hgain  for  simplicity  the  bias,  or  mean  value  of  8e,  will  be 
assumed  to  be  zero  for  the  performance  comparison.  Thus,  again  if  the  SNR  is  rea¬ 
sonable,  the  bearing  error  will  be  small.  Under  this  circumstance  (2-6)  can  be 
approximated  by 


*  S  <  *d  CO®  ©o)  (  ©i  “  ©o> 


This  expression  can  be  cast  into  the  form 


KPF  (6a  -eo) 

where  •  Horton5  phase  factor  t 


(2-9) 


Kpp  s  *  d  cos  6q 


(2-10) 


Equation  (2-9)  shows  that  the  desired  bearing  estimate  can  be  obtained  from  a 
measurement  of  differential  signal  phase.  Putting  it  another  way,  information  con¬ 
cerning  bearing  is  contained  in  signal  phase.  Two  different  ways  of  extracting 
this  bearing  information  are  described  in  the  next  section. 


The  important  fact  in  all  this  is  that  the  phase  estimate  provided  by  the  mea¬ 
surement  of  differential  signal  phase  is  the  resultant  phase  of  the  signal  plus 
accompanying  noise.  This  aspect  is  treated  in  detail  in  subsequent  sections  of 
this  report. 


2-3 


3.0  BASIC  CCHFIGURATIONS 


In  this  section,  the  basic  configurations  of  the  two  bearing  estimators,  whose 
performance  is  being  compared  in  this  report,  are  shown  and  described. 

Figure  3-1  shows  the  closed-loop  bearing  estimator.  Because  of  familiarity, 
it  is  based  on  the  DIFAR-Like  Tracker  (DLT)  analyzed  in  the  DLT  study  report*. 

The  assumed  signal  waveform  and  wavefront  are  as  described  in  the  previous 
section:  a  sinusoidal  signal  and  a  plane  wave.  The  split-array  beamformer  is  also 
described  there. 

To  minimize  the  effect  of  additive  noise  accompanying  the  received  signal, 
bandpass  filters  are  used.  It  is  assumed  in  this  report  that  the  received  signal 
lies  at  the  center  of  the  band  (f q) .  The  quantitative  effect  of  the  bandwidth 
(Af^)  is  determined  analytically  in  this  report. 


A 

0|t0O 


INCOMING 
ACOUSTIC 
PLANE  WAVE; 
ACTUAL  BEARING 

(0|) 


Figure  3-1.  Closed-Loop  Bearing  Estimator 


R2291 


The  phase  detection  function  in  the  DLT  closed-loop  bearing  estimator  is  per¬ 
formed  by  a  clipper  bearing  deviation  indicator  (BDI) .  More  will  be  said  on  this 
subject  in  the  next  section. 

The  remainder  of  the  loop  contains  the  necessary  equaliser  network  for  loop 
stability,  and  the  motor  which  drives  the  steering  comparator  in  the  beamformer. 

The  output  of  the  motor  is  the  bearing  estimate  (0Q) ,  which  is  driven  by  the  sys¬ 
tem  in  such  a  direction  as  to  minimize  the  bearing  error.  The  motor  provides  a 
single  integration,  resulting  in  a  Type-I  loop.  This  is  all  a  description  of  a 
classical  closed-loop  feedback  system.  There  is  no  need  to  dwell  on  these  matters 
in  this  report;  all  analytical  influences  are  already  worked  out  in  the  DLT  study 
report1. 

Figure  3-2  shows  the  basic  configuration  of  the  open- loop  bearing  estimator. 
The  same  assumptions  are  made  covering  the  incoming  signal  and  beamformer,  except 
that  the  main  response  axis  (MRA)  of  each  beam  is  now  fixed.  Probably  a  set  of 
preformed  beams  would  be  used.  Actually,  no  such  system  can  be  truly  open- loop;  at 
least  some  kind  of  loose  control  is  required,  if  only  to  select  the  proper  beam  and 
transfer  between  beams  when  required.  Because  of  the  fixed  beam  MRA,  the  informa¬ 
tion  carried  in  the  differential  phase  is  the  input  bearing  rather  than  bearing 
error.  For  the  purposes  of  this  report,  bearing  rate  of  the  signal  source  is 
assumed  to  be  zero. 

The  configuration  of  Figure  3-2  is  based  on  the  use  of  a  commercial  spectrum 
analyzer:  the  Hewlett-Packard  Model  HP  3582A.  This  piece  of  equipment,  it  turns 
out,  is  inherently  capable  of  performing  a  number  of  the  required  functions  of  the 
bearing  estimator. 

Bandpass  filtering  is  used,  as  it  is  in  the  closed-loop  configuration,  but 
this  time  it  is  provided  by  the  Fast  Fourier  Transformation  (PFT)  contained  in  the 
spectrum  analyzer.  The  parameters  of  a  bin  are,  again,  center  frequency  (fQ)  and 
bandwidth  (Af^j .  The  dependence  of  this  bandwidth  on  FFT  parameters  will  be  dis¬ 
cussed  later.  The  FFT  bin  output  is  conplex;  that  is,  it  appears  as  a  pair  of  com¬ 
ponents:  in-phase  (u)  and  quadrature  (v) .  There  are  actually  two  such  pairs  in 
the  system,  one  for  the  left-hand  channel  and  the  other  for  the  right-hand  channel. 
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Figure  3-2.  Open-Loop  Bearing  Estimator 

The  phase-detection  function  chosen  for  the  open-loop  bearing  estimator  is  the 
argument  of  the  cross-spectral  density  between  the  left-half  and  right-half  FFT  bin 
outputs;  this  is  a  function  contained  within  the  spectrum  analyzer.  The  definition 
of  this  quantity  is  given  in  the  next  section. 

The  phase  estimates  (»|i)  are  delivered  as  a  sequence  of  discrete  samples;  these 
samples  are  averaged  arithmetically  in  the  spectrum  analyzer.  This  time  averaging 
is  discussed  in  a  later  section,  and  compared  with  that  occurring  in  the  closed- 
loop  bearing  estimator. 

Finally,  the  time-averaged  phase  estimate^)  is  divided  by  the  Horton  phase 
factor  to  arrive  at  the  bearing  estimate: 


3-3 


R2291 


4.0  PHASE  DETECTION 

In  the  previous  section,  in  which  the  basic  configuration  of  the  two  bearing 
estimators  is  presented,  the  respective  implementations  of  the  phase-detection 
function  appear  to  be  markedly  different.  Nevertheless,  in  the  present  section 
they  are  seen  to  be  mathematically  identical,  with  one  proviso. 

That  proviso  is  that  the  particular  form  of  the  clipper  BDI  used  in  the 
closed-loop  bearing  estimator  be  one  that  employs  flip-flops  to  achieve  a  sawtooth 
characteristic  spanning  all  four  quadrants;  that  is,  +180°.  Such  a  circuit  can  be 
found  in  a  paper  by  Lindgren  et  al®,  for  example.  When  this  circuit  is  used,  the 
resulting  output  is  equal  to  the  differential  phase  (4<)  between  the  voltages 
appearing  in  the  two  (left-  and  right-half)  channels,  after  bandpass  filtering.  In 
equation  form,  these  two  channel  voltages  can  be  expressed  in  two  ways,  as  follows: 

A  cos  (2TrfQt  +  $L)  +  nL(t)  «  cos  (2Trf0t  +  4>L)  (4-1) 

A  cos  (2rrf gt  +  $R)  +  nR(t)  «  CR  cos  Uir^t  +  <j>R)  (4-2) 

where  the  left  sides  are  in  the  form  of  the  sum  of  a  sinusoidal  signal  with  con¬ 
stant  amplitude  (A)  and  true  phase  ($) ,  plus  Gaussian  noise  with  mean  zero  and 
variance  (N) ,  statistically  independent  from  channel  to  channel.  The  right  sides 
are  in  the  form  of  a  composite  voltage  with  slowly  varying  envelope  (C)  and  phase  (<J>) 
The  output  of  a  +180°  BDI,  in  these  terms,  will  be 

4,  +  180°  *  (4>L  -  4>r  +  1800)  mod  360°  (4-3) 

defined  over 

-1800  <  ^  <  +180°  (4-4) 
The  statistical  properties  of  4*  will  be  discussed  in  the  following  section. 

Attention  is  next  turned  to  the  other  implementation  of  the  phase-detection 
function:  the  argument  of  the  cross-spectral  density.  This  quantity  is  estimated 
in  the  following  fashion:  The  complex  coefficients  in  the  two  (left-  and  right- 
half)  FFT  bin  outputs  associated  with  a  particular  frequency  (fQ)  are  multiplied, 
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taking  the  complex  conjugate  of  one  of  the  outputs  prior  to  multiplication.  The 
cross-spectral  density  is  a  complex  number;  in  polar  form*  its  components  are  a 
magnitude,  which  is  the  commonly  used  output,  and  an  angle,  called  the  argument, 
which  is  used  for  the  present  application. 

Now  the  above  word  procedure  will  be  illustrated  mathematically.  The  FFT  bin 
output  components  are  as  follows: 


u_  =  A  cos  $>  +  n_  _  *  CT  cos  <t>T 

L  L  L  f  X  L  L 

vT  *  A  cos  $  _  +  n.  _  *  C_  cos  <fc  T 

L  L  LrQ  L  L 

uR  -  A  cos  *„  ♦  -  CR  cos  *R 

vR  =  A  cos  *R  *  nR_5  -  CR  cos  *R 

where  the  u  quantities,  and  the  I  subscript,  refer  to  the  in-phase  components; 

while  the  v  quantities,  and  the  Q  subscript,  refer  to  the  quadrature  components. 

The  first  operation  is  the  following  complex  multiplication: 


(UL  +  jVL)  (UR  '  jV 


(4-5)  j 
(4-6) 
(4-7) 
(4-8) 


CLe 


CRe 


-j*R 


CI»CR  e 


j«oL  -<pr) 


(4-9) 


The  argument  (polar  angle)  of  this  complex  expression  is,  as  implemented  in  the 
commercial  instrumentation. 


4j  +  1800  =  (£  -  6  +  180°)  mod  360° 

Li  I\ 


(4-10) 


defined  over 


-180  0  <  4,  <  +180° 


which  is  identical  to  the  +180°  BDI  output,  given  by  (4-3)  and  (4-4). 


(4-11) 
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Thus,  the  general  model  £or  the  phase-detection  function  can  be  illustrated  as 
in  Figure  4-1.  The  input  is  (jr)  the  differential  phase  of  the  signal,  while  the 
output  is  (4«)  the  differential  phase  of  the  signal  plus  noise. 


*  +  180*  - 
($L  -  +  180*) 

MOD  360* 


PHASE 

DETECTOR 


«l»  +  180*  - 
(4>l  -  4>  r  +  1 80*) 
MOD  360* 


Figure  4-1.  Definition  of  Phase  Quantities 
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5.0  PHASE  STATISTICS 

In  the  previous  section,  the  outputs  of  the  phase  detectors  associated  with 
the  closed-  and  open-loop  bearing  estimators  are  shown  to  be  mathematically 
identical.  In  the  present  section,  the  statistics  of  this  common  output  are 
shown.  The  history  and  results  of  the  analysis  developing  these  statistics  are 
presented  here;  the  details  are  given  in  Appendix  A. 

An  earlier  analysis  appears  in  the  DLT  study  report1-.  There  the  configura¬ 
tion  of  the  well-known  clipper  BDI,  which  operates  in  the  first  and  fourth  quad¬ 
rants  of  phase — that  is,  +90° — is  shown.  A  rather  detailed  analysis  is  performed 
in  that  report  to  obtain  the  statistics  of  the  phase  estimate  provided  by  this  par¬ 
ticular  circuit.  This  analysis  is  based  on  the  probability  density  function  (pdf) 
of  the  phase  of  a  sine  wave  plus  Gaussian  noise,  derived  by  Middleton  in  19487. 

Q 

Later”,  he  presented  a  Fourier  series  representation  of  the  same  pdf.  Using  this 
form,  and  assuming  the  additive  noise  terms  in  the  two  input  channels  to  be  statis¬ 
tically  independent  of  one  another,  the  analysis  in  the  DLT  study  derives  Fourier 
series  expressions  for  the  moments  of  the  differential  phase  estimate  as  a  function 
of  the  input  signal-to-noise  ratio  (SNR)  after  bandpass  filtering. 

Now,  it  was  also  required  in  the  prosecution  of  the  DLT  study  to  determine  the 
effect  on  systems  performance  of  employing  the  other  form  of  clipper  BDI,  extended 
to  cover  all  four  quandrants  (+180°).  This  is  the  form  considered  in  the  present 
report.  The  resulting  series  are  presented  in  the  DLT  study  report,  but  the  deri¬ 
vation,  which  is  parallel  to  that  presented  for  the  conventional  (+90°)  BDI,  is  not 
shown  there. 

Another  feature  of  the  phase  statistics  derived  in  the  DLT  study  is  their 
form.  Although  the  mean  and  variance  of  the  differential  phase  estimate  («^)  are 
both  functions  of  the  true  differential  phase  (>&) ,  as  well  as  the  input  SNR,  the 
specific  statistics  chosen  in  that  study  are  the  slope  of  the  mean  curve  of  4,  ver¬ 
sus  'l'  evaluated  at  the  origin: 

9*  'I'—O 
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and  the  standard  deviation  (square  root  o£  the  variance)  of  41,  also  evaluated  at  ♦ 
equal  to  zero. 

Thus,  in  Appendix  A  of  the  present  report,  two  additional  actions  are  taken: 

a)  The  derivation  of  the  phase  statistics  for  the  extended  (+180°)  BDI  is 
shown  explicitly.  The  work  is  similar  to  that  by  others9-1*.  Because  of 
the  mathematical  identity  of  the  outputs,  these  statistics  apply  also  to 
the  argument  of  the  cross-spectral  density. 


b)  The  selected  statistics  are  now  the  mean  and  standard  deviation  of  the 
phase  estimate  (41)  versus  the  actual  differential  phase  C&)  for  various 
values  of  SNR. 

The  results  of  the  derivation  in  Appendix  A  are  as  follows  (the  units  of  v 
and  ^  are  in  radians  throughout):  The  first  moment,  or  mean,  of  the  phase  estimate 
is  given  by 


uu 

*)  "  2  X 


+  1 


m  *=  1 


B  (z)  sin  m* 
in 


The  second  moment,  or  mean  square,  of  the  phase  estimate  is  given  by 


«l»  2(* 


-r~<  2,  ^ 


m  +  1 


m  =  1  m 


- -  B  (z)  cos  m 

2  m 


where  ■  coefficient,  expressible  in  various  alternative  ways,  as  shown  in 
Appendix  A. 

The  quantity  z  is  defined  as 


2  z 

*o  "  T 

where  a*  ■  Middleton's8  symbol  for  input  signal-to-noise  power  ratio  (SNR) 


In  terms  of  system  parameters, 


5-1 


R2291 


a 


2 

0 


*0**1 


(5-4) 


where 


A  *  amplitude  of  sinusoidal  signal 

Nq  *  spectral  density  of  noise  (assumed  flat  over  the  pass  band) 

Af^  *  effective  noise  bandwidth  of  input  bandpass  filters. 

During  the  course  of  the  study,  it  was  found  that  it  is  quite  possible  to  pro¬ 
gram  the  numerical  calculation  of  (5-1)  and  (5-2) ,  together  with  the  coefficients 
on  a  programmable  hand  calculator:  specifically  the  Hewlett-Packard  HP-41C  calcu¬ 
lator.  Three  memory  modules  are  required,  and  the  assignment  of  150  data  registers 
is  sufficient.  The  calculation  is  relatively  slow,  but  definitely  feasible.  The 
details  of  the  program  are  being  submitted  to  the  Government  under  separate  cover. 


Figure  5-1  is  a  plot  of  the  results  of  the  calculation.  The  units  of  the  plot 
are  in  electrical  degrees  rather  than  radians.  Shown  are  the  mean  and  standard 
deviation  (called  "SIGMA*  in  the  figure)  of  the  phase  estimate  (4/)  versus  the  act¬ 
ual  differential  phase  (30  for  several  values  of  SNR  in  dB. 

The  features  to  be  noted  are  as  follows:  The  moments  of  4/  are  expressed  in 
equations  (5-1)  and  (5-2)  in  Fourier  series  form,  hence  are  periodic  in  3>.  These 
two  equations  are  derived  in  Appendix  A  based  on  a  sawtooth  characteristic  spanning] 
+180°:  thus  it  can  be  seen  that  the  period  of  the  moments  is  360°.  The  mean,  given 
by  (5-1),  appears  as  a  summation  of  sine  functions,  hence  is  an  odd  function  of  3'. 
On  the  other  hand,  the  standard  deviation  is  an  even  function  of  3',  for  the  follow¬ 
ing  reason:  In  equation  form 


V*|*>  * 

The  first  term  under  the  radical  is  given  by  (5-2) ,  which  appears  as  a  summation  a 
cosine  functions,  hence  is  an  even  function  of  3».  The  second  term  is  the  square  o 
an  odd  function,  therefore  ar.  even  function.  The  difference  is  even,  and  so  is  th 
positive  square  root  of  the  difference. 


J*2  (*|z) 


W  z) 


(5-1 


54 


80' 


Figure  5-1.  Moments  of  Phase  Estimate,  ijj 


The  mean  curve  has  a  central  linear  portion,  whose  extent  is  greater  when  the 
SNR  is  higher.  The  slope  of  this  linear  portion,  which  is  important  in  what  fol¬ 
lows,  is  a  function  of  SNR,  becoming  smaller  as  SNR  decreases.  The  bias  is  the 
difference  between  the  mean  of  and  the  true  valued;  this  is  shown  in  the  fig¬ 
ure.  Its  value  is  always  zero  or  negative,  because  4>  is  always  less  than  or  equal 
to 


The  standard  deviation  curve  has  a  flat  central  region,  whose  extent  again  is 
greater  when  the  SNR  is  higher.  The  value  of  the  standard  deviation  in  this 
region,  which  is  also  important  in  what  follows,  is  a  function  of  SNR,  becoming 
larger  as  SNR  decreases. 

The  practical  system  effects  of  these  attributes  of  the  phase  estimation  sta¬ 
tistics  become  more  clear  in  following  sections  of  this  report. 
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6.0  LOCATION  AMD  NATURE  OF  ERROR 


Although  the  two  bearing  estimators  under  study  in  this  paper  possess  phase 
detection  implementations  with  the  same  statistics  (which  are  presented  in  the  pre¬ 
vious  section  of  this  report) ,  there  is  still  a  fundamental  theoretical  difference 
between  the  estimators:  the  location  and  nature  of  the  error  used  as  a  (negative) 
measure  of  performance. 


In  the  closed- loop  bearing  estimator,  the  bearing  error,  which  is  the  differ¬ 
ence  between  the  output  bearing  estimate  (0Q)  and  the  input  true  bearing  (0^) , 
can  be  transformed  by  the  Horton  phase  factor  of  the  beamformer  into  a  correspond¬ 
ing  phase  error  at  the  input  to  the  phase-detection  function;  that  is,  an  error  in  ^ 
In  this  report,  this  normalization  is  in  fact  made,  to  eliminate  the  necessity  of 
specifying  operating  frequency  and  aperture  size.  However,  the  source  of  the  error 
is  to  be  found  in  the  output  statistics  of  the  phase  detection  function;  that  is, 
the  statistics  of  4/ . 


Now,  to  determine  the  degree  of  error  in  S'  caused  by  the  statistics  of  4, 
requires  a  loop  analysis;  strictly  speaking,  one  that  takes  into  account  the  non¬ 
linear  nature  of  the  phase  detection  function.  Section  4  of  the  dlt  study  report1 
presents  the  history  and  present  state  of  the  art  of  nonlinear  loop  analysis.  As 
pointed  out  there,  the  use  of  an  analytical  technique  involving  the  Pocker-Planck 
equation  leads  to  a  probability  density  function  (pdf) ,  thence  to  the  desired 

moments,  in  a  particularly  sinple  loop  configuration  (sinpler  than  the  DLT  loop 

12 

under  consideration  here) ,  the  result  is  the  Tikhonov  distribution: 


pW 


a  cos  St 
e 

2  ttI Q(a) 


(6-1) 


where 


a  *  the  signal-to-noise  ratio  in  the  bandwidth  of  the  loop, 

Iq  *  modified  Bessel  function  of  the  first  kind,  zero^  order. 
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This  Tikhonov  distribution  has  a  relatively  simple  appearance,  tempting  one  to 
make  use  of  it,  despite  its  lack  of  strict  applicability  to  the  DLT  loop,  as  an 
approximation  to  account  for  the  large-error  domain.  In  the  end,  however,  the 
decision  has  been  made  to  resort  once  again  to  the  linearized  loop  analysis  used  in 
the  DLT  study.  (In  making  such  a  decision,  it  is  necessary  to  be  aware  at  the  out¬ 
set  that  a  linearized  analysis  does  not  speak  to  the  problems  of  acquisition  and 
cycle  slipping.) 


In  the  linearized  analysis,  the  slope  of  the  central  linear  portion  of  the 
mean  of  ^  is  regarded  as  the  gain  factor  of  the  phase  detector.  The  symbol  (KpD) 
will  be  used  in  this  report,  although  the  quantity  so  defined  differs  from  that 
with  the  same  symbol  in  the  DLT  study  report  by  a  factor  of  tt  for  the  +180°  BDI. 
The  definition  in  the  present  report  is 


K 


PD 


(6-2) 


When  this  is  done,  the  normal  action  of  the  closed  feedback  loop  is  to  yield 
the  following  equation1: 


<V°IZ)  /Afe 
kpd  J  *V2 


(6-3) 


where 


Af@  »  effective  noise  bandwidth  of  the  closed  loop, 

Af^  m  effective  noise  bandwidth  of  the  input  bandpass  filters. 

The  other  factors  have  already  been  defined  in  this  report. 

At  this  point,  it  is  well  to  mention  that,  in  the  absence  of  a  bearing  rate, 
the  closed-loop  bearing  estimator  produces  an  unbiased  estimate.  Put  another  way, 
this  Type-I  loop  (single  integration)  produces  no  lag  error  in  the  absence  of  a 
bearing  rate. 
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Turning  to  the  open-loop  bearing  estimator,  the  output  bearing  estimate  (6o) 


is  given  by 


where 


(  )=  arithmetic  average  over  K  samples, 

KpF  *  Horton  phase  factor,  defined  earlier. 

As  before,  normalization  is  applied  by  removing  the  Horton  phase  factor,  so 
that  operating  frequency  and  aperture  size  need  not  be  specified.  The  quantity  to 
be  studied  becomes  simply  (^),  the  time-averaged  phase  at  the  output  of  the  phase 
detector.  The  phase  estimate  at  the  output  of  the  phase  detector  is  in  general 
biased,  in  the  previous  section,  the  mean  (i£)  is  plotted  versus  the  true  phase  (>£) 
for  various  values  of  SNR.  It  can  be  seen  that  there  is  a  difference  between  ^ 
and  ^  (that  is,  a  bias)  for  nonzero  'i  and  finite  SNR.  The  time-averaging  process 
will  reduce  the  standard  deviation  of  the  estimate,  but  will  not  affect  the  bias. 

For  the  purpose  of  the  performance  comparison  in  this  report,  the  question  of 
bias  will  be  ignored.  It  is  assumed  that  the  sound  source  is  on  the  beam  main 
response  axis  (MRA) ,  so  that  the  actual  differential  phase  ('£)  is  zero.  Figure  5-1 
shows  that  this  results  in  the  mean  phase  estimate  (i|>)  going  to  zero,  independent 
of  SNR;  hence  the  bias,  which  is  the  difference  between  ^  and  ty,  is  also  zero. 

Having  assumed  away  the  bias,  the  expression  for  error  is  in  terms  of  the 
standard  deviation: 


Equation  (6-3)  predicts  the  theoretical  performance  of  the  closed-loop  bearing 
estimator,  while  equation  (6-5)  does  the  same  for  the  open-loop  bearing  estimator. 
Now,  these  two  equations  are  difficult  to  interpret  and  compare  in  their  present 
forms.  Therefore,  they  will  be  manipulated  further  in  the  next  section. 
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7.0  COMPARATIVE  PERFORMANCE 


In  this  section,  the  £inai  performance  expression  will  be  developed  and  the 
results  of  a  numerical  evaluation  of  these  expressions  will  be  plotted  to  show  the 
comparative  performance  of  the  closed-  and  open-loop  bearing  estimators. 

Both  preliminary  performance  equations,  (6-3)  and  (6-5),  include  the  idea  of 
the  reduction  of  fluctuations  by  time-averaging  over  a  fixed  time  interval  called 
the  observation  tin®,  which  can  be  called  T Qbs .  The  resulting  standard  deviation 
in  each  case  is  inversely  proportional  to  the  square  root  of  T  For  further 

universality  of  the  results,  it  is  possible  to  eliminate  the  necessity  of  specify¬ 
ing  a  value  for  this  observation  time  by  evaluating  the  product  of  the  standard 
deviation  and  the  square  root  of  TQbs.  This  idea  is  not  original  with  this 
report;  Fischer13  employed  it  in  an  earlier  comparative  study. 


The  closed-loop  bearing  estimator  is  taken  up  first.  Multiplying  (6-3)  by  the 
square  root  of  the  observation  time  yields  the  following: 


<T 


* 


obs 


(7-1) 


The  following  rule  of  thumb  is  suggested,  relating  the  equivalent  observation 
time  to  the  effective  noise  bandwidth  of  the  analog  loop: 


During  the  course  of  this  study,  the  validity  of  (7-2)  has  been  tested  empirically 
using  real  data  from  the  DLT  Engineering  System  Test  Model  (ESTO) . 

On  the  theoretical  side,  the  effective  noise  bandwidth  is  computed  for  the 
specific  system  used  for  taking  the  data  in  question  by  making  use  of  the 
equation1 


Af 


e 


4T 


0.3 
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Substituting  (7-2)  into  (7-1)  yields 


<r 


Mr 


(7-4) 


Next,  attention  is  turned  to  the  open-loop  bearing  estimator.  Equation  (6-5), 
repeated  here  for  convenience. 


°V°  |*> 

7^ 


(6-5) 


is  written  in  terms  of  the  quantity  K,  the  number  of  samples  of  ^  entering  the  time 
averaging  process.  There  is  implicit  in  (6-5)  the  assumption  of  statistically 
independent  samples.  This  matter  deserves  further  discussion. 


In  the  first  place,  the  samples  in  question  are  taken  from  an  FFT  bin,  the 
input  to  which  is  a  sequence  of  input  data  of  time  duration  Tfffc.  Now,  a  process 
called  "overlapping"  is  sometimes  used,  wherein  the  first,  say  fifty  percent,  of 
the  input  sequence  associated  with  one  output  sample  is  identical  to  the  last  fifty 
percent  of  the  input  sequence  associated  with  the  previous  output  sample.  In  such 
a  case,  the  output  samples  would  be  correlated  with  each  other.  Such  overlapping 
is  not  intended  in  the  open-loop  bearing  estimator,  and  is  not  provided  in  the 
HP-3582A  spectrum  analyzer. 

In  the  absence  of  overlapping,  the  output  samples  tend  to  be  uncorrelated  with 
each  other  to  a  large  degree,  particularly  if  a  large  number  of  input  samples  (for 
example,  1024)  is  used  in  the  input  sequence.  There  may  be  a  small  "edge  effect* 
representing  the  correlation  between  input  samples  at  the  beginning  of  one  sequence 
and  input  samples  at  the  end  of  the  previous  sequence.  This  effect  should  be  quite 
negligible  in  practice. 

Thus  the  expression  for  the  number  K  is  simply 


K 


(7-5) 
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Therefore 

"<+> 

In  order  to  have  a  common  basis  for  performance  comparison,  it  would  be  desir¬ 
able  to  write  (7-6)  in  terms  of  the  input  bandwidth, Af  ^  This  bandwidth  has 
been  specified  as  the  effective  noise  bandwidth  of  the  input  bandpass  filter,  in 
this  case  the  FFT  bin. 

Harris14  has  published  a  comprehensive  treatment  of  the  quantitative  effect 
of  "window"  shape  on  a  number  of  FFT  parameters,  including  effective  noise  band¬ 
width.  A  window  is  a  time-domain  weighting  function  used  to  improve  the  qualities 
of  the  spectral  estimates.  The  effective  noise  bandwidth  can  be  expressed  as  fol¬ 
lows: 


*  ,  (0  z) 


(7-6) 
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All  that  remains  for  a  performance  comparison  is  a  numerical  evaluation  of 
(7-4)  for  the  closed-loop  bearing  estimator  and  (7-9)  for  the  open- loop  bearing 
estimator.  The  right-hand  sides  of  the  two  equations  are  identical,  except  for  the 
phase  detector  gain  factor  (Kp0) .  They  both  require  the  standard  deviation  of 
the  phase  estimate  given  by  (5-6)  evaluated  where  the  true  differential  phase  MO 
is  zero.  In  addition,  the  quantity  (KpD) ,  given  by  (6-2) ,  is  expressed  in  series 
form  and  evaluated  for  use  in  (7-4) .  The  quantity  is  also  evaluated  in  the  OLT 
study  report^-,  but  differs  from  the  corresponding  quantity  in  the  present  report 
by  a  factor  of  it  ,  because  of  definition. 

Figure  7-2  is  a  plot  of  the  results  of  the  above-described  calculation.  The 
appropriate  normalized  phase  error  used  as  a  (negative)  measure  of  performance  in 
this  report  is  plotted  versus  signal-to-noise  ratio  in  the  half  beam  and  in  a  1-Hz 
band,  this  latter  for  further  universality  of  the  results.  In  decibels,  the 
abscissa  is 

10  log-^ji- 

The  input  bandwidth  (M^)  appears  as  a  parameter  of  these  curves. 

By  inspection  of  Figure  7-2,  it  can  be  seen  that  the  performance  of  the  two 
bearing  estimators  is  different  at  low  signal-to-noise  ratios,  because  of  the  pres¬ 
ence  of  the  phase  detector  gain  factor  in  the  closed-loop  expression  (7-4) .  This 
factor  causes  the  phase  error  to  increase  as  the  signal-to-noise  ratio  decreases. 
Physically,  this  increase  can  go  on  only  so  far  before  track  is  broken  and  cycle¬ 
slipping  occurs,  a  phenomenon  not  considered  explicitly  in  this  report.  As  men¬ 
tioned  earlier,  a  different  and  more  involved  analytical  procedure  is  required  to 
account  for  this  phenomenon.  If  such  a  procedure  were  to  be  used,  probably  the 
curves  would  be  more  similar  to  each  other  in  this  low  SNR  region. 

The  phase  error  for  the  open-loop  estimator,  on  the  other  hand,  can  be  seen  to 
approach  a  constant  asymptote,  whose  value  depends  on  the  input  bandwidth  , 

as  the  SNR  decreases.  The  asymptotic  value  is  that  associated  with  the  uniform 
phase  distribution  characteristic  of  pure  Gaussian  noise. 
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Figure  7-2.  Comparative  Performance 


The  important  conclusion  of  the  paper  is  arrived  at  by  an  examination  of  the 
performance  curves  at  high  signal-to-noise  ratios.  It  can  be  seen  that  all  curves 
approach  each  other,  and  the  performance,  defined  as  it  is  in  this  paper,  becomes 
independent  of  input  bandwidth,  but,  more  importantly,  becomes  the  same  for  the 
closed-loop  and  open-loop  bearing  estimators. 

Thus,  the  decision  as  to  which  of  the  two  estimator  concepts  to  use  should  be 
based  on  practical  implementation  considerations  rather  than  theoretical  perfor¬ 
mance.  Operational  system  studies  will,  of  course,  need  to  address  the  question  of 
source  motion  and  statistical  estimate  bias  deliberately  excluded  from  considera¬ 
tion  in  this  report. 
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APPENDIX  A 

DERIVATION  OP  PHASE  STATISTICS 


In  this  appendix,  the  statistics  of  the  estimate  of  differential  phase  pro¬ 
vided  by  the  phase-detection  function  of  the  two  bearing  estimators  of  this  report 
will  be  derived.  In  the  main  text  it  is  shown  that  the  two  phase-detection  imple¬ 
mentations,  while  outwardly  quite  different,  provide  mathematically  identical  phase 
estimates.  Also,  it  turns  out  that  much  of  the  derivation  has  already  been  given 
in  Section  6.8  of  the  DLT  study  report^-.  In  addition,  others®-^--*-  have  pub¬ 
lished  similar  results. 

The  specific  quantities  derived  here  are: 

a)  Mean  of  the  differential  phase  estimate  as  a  function  of  the  true  differen¬ 
tial  phase,  and  input  signal-to-noise  ratio. 

b)  Standard  deviation  of  the  differential  phase  estimate  as  a  function  of  the 
true  differential  phase,  and  input  signal-to-noise  ratio. 

c)  First  derivative  (slope)  of  the  phase  estimate  mean  with  respect  to  the 
true  differential  phase,  evaluated  where  the  latter  independent  variable  is 
zero,  versus  input  signal-to-noise  ratio.  This  quantity  is  called  the 
phase-detector  gain  factor  in  this  report  (with  a  caveat  about  its  scale 
factor) . 

The  derivation  starts  with  the  expression  for  the  vth  moment  of  the  differ¬ 
ential  phase  estimate: 

—  r+7T  r+7r 

(WZ)  =  J  J  *V<*L  -  V  P(*L  |VZ)  •  p(4>r  I  Vz)  d\ 

—  IT  —  IT 

(A-l) 


A-l 


R2291 


where 

4<  =  differential  phase  estimate 

<t>LfR  ■  phase  of  the  signal-plus-noise  voltage  in  the  (left,  right) 
channel 

$L,R  *  phase  of  the  signal  in  the  absence  of  noise  (true  phase)  in  the 
(left,  right)  channel 

p(4>|$;z)  ■  conditional  probability  density  function  (pdf)  of  4> ,  given  values 
of  <£  and  z  (defined  below) 

z  *  input  signal-to-noise  ratio  (not  in  dB) ,  divided  by  2. 

At  this  point  it  is  well  to  point  out  that  the  form  of  (A-l)  with  the  prodjct 
of  two  pdf  factors  appearing,  implies  that  the  noises  in  the  two  (left  and  right) 
channels  are  statistically  independent  of  each  other.  This  is  a  common  assumption 
in  studies  of  this  type,  although  it  is  recommended  that  this  assumption  be  relaxed 
in  any  follow-on  studies. 

There  is  more  to  the  story  than  independence,  however.  The  fact  that  the  two 
pdf  factors  appear  without  subscripts  implies  that  they  are  functionally  identi¬ 
cal.  Also  the  same  value  of  the  parameter  z  appears  in  both  factors,  indicating  an 
identical  input  signal-to-noise  ratio  (SNR)  existing  in  both  channels.  The  DLT 
study  report*  went  into  the  question  of  unequal  SNR  values,  but  the  present  study 
does  not. 

The  next  step  is  to  specify  the  functional  form  of  ^ .  In  the  main  text,  the 
following  is  established  for  both  phase  detector  implementations: 

4»  +  180°  -  (4>l  -  4>r  +  180')  mod  360'  (A-2) 

This  rather  compact  algebraic  expression  describes  the  periodic  sawtooth  character¬ 
istic  pictured  in  Figure  A-l.  This  figure  resembles  Figure  6-4b  of  the  DLT  study 
report*,  but  the  vertical  axis  has  a  different  scale  factor. 


SO' 


(*L“  *R> 

Figure  A-l.  Sawtooth  Characteristic 


Equation  (A-2)  cannot  be  used  in  (A-l)  in  its  present  form.  Since  the  proba 
bility  density  function  (p)  will  be  seen  to  be  expressible  in  Fourier  series  form 
it  will  be  mathematically  convenient  to  express  41  in  the  same  way.  Because  the 
function  is  periodic,  this  form  is  quite  appropriate  and  natural. 


Formula  (509)  of  Jolley15  furnishes  the  required  series,  after 
tary  manipulation: 


some  elemen- 


+.2  5  -4 

m  ■  1 


*  ♦  1 


sin  *(<j>L  -  <t>R) 


It  will  be  noted  that  this  series  is  not  the  same  as  (6-58)  of  the  DL T  study 
report1.  The  latter  series  applies  to  the  triangular,  or  +90* ,  characteristic  of 
Figure  6-4a  instead. 
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Now,  a  secies  foe  4*  is  not  enough;  one  is  required  for  4^  as  well.  Again, 
Jolley  furnishes  the  required  series,  this  time  in  Formula  (517),  after  some  ele¬ 
mentary  manipulation: 


-  00 

4-*  2* 


m  +  1 


m  ■  1  m 


cos  m(4>L  -  4>R) 


(A-4) 


Now,  (A-3)  and  (A-4)  are  written  in  terms  of  the  difference  between  <|>L  and 

4>r  separately.  It  will  therefore  be  convenient  to  rewrite  (A-3)  and  (A-4)  in  a 
similar  fashion,  through  the  use  of  trigonometric  identities: 


r  (.ii"1  +  1 

4i  «  2  2j  -  sin  COS  m<^R 


+  1 


m  -  1 


cos  m4>  sin  m<t> 

L  K 


,  2  ^  m  +  1 

■  — - 4  2^  *  -  cos  m4>L  cos  m<(>R 

tt  *  1  m 

+  1 

+  4  2^  ■^^2 -  8in  m<^L  8in  m<^R  (A-6) 

m  »  1  nr 

It  is  now  possible  to  substitute  (A-5)  and  (A-6)  into  (A-l) ,  interchanging  the 
summations  and  integrations,  which  is  legitimate  here: 


2  I  ^ 

m  «  1 


m  +  1 


sin  m$L  cos  m4>R 


m  +  1 


^  — “ —  cos  m4»L  sin  m^ 


m  »  1 


— r  2  r’  .  +  1 _ 

^  m - 4  2^  '  -  COS  008 

3  m  «  1  mz 


V  (_i\ro  +  1  _ 

+  4  }  >  -  sin  m4»  sin  «4 

_4—i  ,  _2 

m  ■  1  m 


(A-8) 
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where 

_  5  r+" 

sin  n>4>  J  pM>|$;z)  sin  m4>  d4>  (A-9) 

—  IT 

_  =  ^  +  TT 

cos  m<t>  /  p(4>|4>;z)  cos  ra4>  d4>  (A-10) 

-IT 

The  specific  expression  for  (p) ,  the  probability  density  function  of  the  phase 
of  a  sinusoid  plus  Gaussian  noise,  is  available,  having  been  derived  by 
Middleton''0.  It  is  not  necessary,  however,  to  write  the  expression  explicitly, 
since  the  results  of  the  integrations  in  (A-9)  and  (A-10)  also  have  been  given  by 
Middleton8 : 
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.m  +  1 


ST1  /—i  \  A  o 

vp  »  2  2_  - - B^  (z)  8in  ®*fc,  008  m\ 


m  ■  1 


y  L-ij_ 


m  +  1 


m  «  1 


B  (z)  cos  m$r  sin  m$_ 
nt  l  k 


(A-14) 


2  ir^  V  (-d®  +  ^  2 

d)  *  — r - 4  /  J — ‘ — - -  B  (z)  cos  m$  cos  m$D 

3  _  ^  .  _2  TO  li  K 

in  *  l  m 


uu  , 

ir  .  ,.m  +  1  2 

/  *  ■  *  - -  B  z(z)  sin  m$>  sin  m<fc_ 

i— i  ,  2  III  la  R 

m  ■  1  m 


(A- 15) 


Using  the  trigonometric  identities  in  reverse. 


tu 

4>  2  ^  sin 

m  1 


m  +  1 


(A-16) 


2  ir2  V  f-l)m  +  1  2 

4<  *  — r - 4  £  J ^ -  Bm  (z)  cos  m^ 

m  *  1  mz 


(A- 17) 


where 


^  *  -  <J> 

L  R 


(A- 18) 


The  following  quantities  have  been  calculated  numerically  in  this  report: 


iM’S'lz) 


or.  (’J'lz) 
4- 


,  2  _  -2 

vu  4* 


Kpn  (*) 


This  last  quantity,  called  the  phase  detector  gain,  differs  from  the  quantity 
by  the  same  name  in  the  DLT  study  report  by  a  scale  factor  of  v.  in  the  present 
report,  the  specific  expression  is 
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00 

KpD  -  2  X  (-I)”  +  1  Bm2  <*>  (A-19) 

m  ■  1 

Earlier  in  this  appendix,  it  was  mentioned  that  alternative  forms  exist  for 
the  expression  (A-13)  for  the  coefficient  Bm.  The  confluent  hypergeometric  func¬ 
tion  can  be  expressed  in  terms  of  the  Bessel  function,  with  the  following  result: 

(z)  +  1  „  ,  ,  (z)l  (A- 20) 

HI  T  1 
2 

where  In  (  )  «  modified  Bessel  function  of  order  n. 

For  numerical  calculations,  there  exist  convenient  recursive  relationships  for 
calculating  all  the  Bessel  functions,  hence  all  the  coefficients  (B^)  in  terms  of 
the  following  basic  functions: 

1q  and  1^ 

I1/2  and  *3/2 

The  values  for  (IQ)  and  (1^)  can  be  calculated  using  equations  9.8.1  through 
9.8.4  of  the  National  Bureau  of  Standards  handbook16.  There  is  a  sufficient  num¬ 
ber  of  registers  in  the  HP-41C  calculator,  with  three  memory  modules  installed,  to 
have  the  necessary  coefficients  for  the  NBS  equations  prestored. 

The  other,  fractional  order,  modified  Bessel  functions,  (Ijy2)  and  <IV2^ 
are  expressible  in  terms  of  the  simple  exponential  function,  which  is  provided  in 
the  calculator. 

As  already  mentioned,  the  details  of  the  program  to  carry  out  the  recursive 
relationships  to  evaluate  the  coefficients,  and  the  subsequent  summation  of  the 
series,  are  being  provided  to  the  Government  under  separate  cover. 
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APPENDIX  B 

DIFFERENTIAL  PHASE  vs  INDIVIDUAL  PHASE 


Early  in  this  study,  it  was  recognized  that  it  is  necessary  to  have  phase  sta¬ 
tistics  in  order  to  predict  system  performance.  Prior  to  the  theoretical  develop¬ 
ment  of  Appendix  A  and  the  subsequent  numerical  calculation,  the  idea  was  conceived 
of  obtaining  phase  statistics  by  a  Monte  Carlo  computer  simulation. 

Such  a  simulation  was  in  fact  performed,  initially  by  NUSC  New  London  Labora¬ 
tory,  and  then  completed  by  Sonalysts,  Inc17.  The  specific  task  was  aimed  at 
comparing  three  different  averaging  techniques: 

I.  <<*>>  -  tan"1 

II.  <*>  -  <tan_1  -JJ-) 

III.  <*)-  tan'1  <-^-> 

where  the  corner  brackets  signify  a  simple  arithmetic  average  over  a  finite  number 
(K)  of  independent  samples.  Several  statements  can  be  made  about  this  study  of 
time-averaging  techniques: 

a)  The  full  details  appear  in  the  Sonalysts  report17. 

b)  In  that  report,  it  is  found  that  Technique  III  is  inadequate.  The  reason¬ 
ing  is  given  there. 

c)  For  the  same  amount  of  averaging  (number  of  independent  samples) ,  Tech¬ 
niques  I  and  II  are  comparable  at  high  signal-to-noise  ratios  (SNR) ,  but 
Technique  II  is  superior  at  low  SNR  where  the  phase  distributions  approach 
the  uniform  distribution  characteristic  of  Gaussian  noise  alone. 

d)  In  fact,  the  selected  bearing  estimator  configurations  utilize  Technique 
II.  For  the  HP3582A  spectrum  analyzer  see  Reference  18. 
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Now,  having  performed  the  Monte  Carlo  simulation,  it  seemed  that  the  results 
could  then  be  used  for  the  comparative  performance  predictions  for  the  open-  and 
closed-loop  bearing  estimators. 

The  first  fact  to  be  noted  is  that  the  Monte  Carlo  simulation  was  of  the  phase 
of  a  sinusoid  plus  Gaussian  noise;  in  other  words,  the  phase  in  one  channel.  The 
phase  detectors  in  the  bearing  estimators  provide  an  estimate  of  differential  phase 
between  two  channels.  It  is  thus  necessary,  if  the  Monte  Carlo  simulation  results 
are  to  be  used  for  performance  prediction,  to  relate  the  statistics  of  differential 
phase  to  those  of  individual  phase. 

The  initial  (mistaken)  theory  was  as  follows:  The  differential  phase  estimate 
provided  by  the  two  phase  detectors  was  assumed  to  be 

4*  -  6l  -  4>R,  (B-1) 

a  simple  subtraction.  If  (B-1)  were  true,  the  mean  and  variance  would  be  as  fol¬ 
lows: 


*  4>($L  J  SNR)  -  4>  (4R;  SNR) 

(B-2) 

%2  =a4,2  ($L;  SNR)  +0^2  ($R;  SNR) 

(B-3) 

$L  *  **  $R  ** 

*  *  ~  *R  *  24 

4^-24  (4/2;  SNR) 

(B-4) 

o-^2  -  20-4,2  (*/2>  SNR) 

(B-5) 
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Equations  (B-4)  and  (B-5)  would  provide  the  required  relationship  between  the 
statistics  of  the  differential  phase  on  the  one  hand  and  individual  phase  on  the 
other,  if  (B-l)  were  true.  Unfortunately,  however,  (B-l)  is  not  true. 

It  has  been  established  in  the  main  text  that  the  correct  expression  for  the 
differential  phase  estimate  provided  by  the  phase  detectors  in  the  two  bearing 
estimators  is  as  follows: 

^  +  180°  =  (4>l  -  4>r  +  180°)  mod  360°  (B-6) 

where  the  180°  term  centers  the  domain  of  definition  about  zero,  and  the  mod  360° 
operation  causes  the  domain  to  extend  over  all  four  quadrants  (360°): 

-180° <  ^  s  +  180° 

The  difference  between  the  correct  expression  (B-6)  and  the  simpler  expression 
(B-l)  resides  in  a  phenomenon  called  "phase  wraparound". 

The  correct  differential  phase  statistics  are  those  derived  in  Appendix  A. 

The  statistics  are  related  to  those  of  individual  phase  in  a  simple  fashion  only 
when  true  phase  is  near  zero  and  the  signal-to-noise  ratio  is  high;  these  condi¬ 
tions  ensuring  a  low  probability  of  "wraparound". 
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